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Abstract 



The propagation of extremely energetic v^s through the Earth is studied by 
means of a detailed Monte Carlo simulation. All major mechanisms of v T interactions 
and r energy loss as well as all its relevant decay modes are properly taken into ac- 
count. The detection possibility of v T J s emerging from Earth in forthcoming neutrino 
telescopes is addressed. 
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1 Introduction 



Cosmic neutrinos with energy > 10 12 eV are a new and exciting observational window 
for cosmic ray physics and astrophysics. Indeed, many interesting models predict the ex- 
istence of cosmic neutrinos with an energy up to 10 18 — 10 20 eV (extreme energy), where 
charged particles are expected to be cut-off by the so called GZK mechanism Q. The 
aforementioned threshold of 10 12 eV is necessary for a detection of diffuse flux, because 
of the large atmospheric neutrino background. 

The detection of such extreme energy neutrinos is challenging as it demands a very large 
sensitive area. The experiments of the upcoming generation plan to use underwater- 
underice Cerenkov detectors || and large field-of-view atmospheric fluorescence detec- 
tors ||. One of the best evidence of cosmic neutrinos would be the detection of upstream 
showers/particles emerging from the Earth. For this kind of events, the atmospheric 
muon and primary charged cosmic ray background would be completely suppressed. On 
the other hand, this signature can hardly be observed at extreme energy because the rise 
of weak cross sections entails the opacity of the Earth with respect to neutrino propaga- 
tionl,|. 

The problem of neutrino propagation through the Earth can be understood on the basis 
of neutrino interactions with matter. For neutrino energy well above few GeV, interac- 
tions with conventional matter are described by neutrino-nucleon deep inelastic scattering 
and neutrino-electron scattering ||. The former can be classified in two main categories: 
charged current (CC) v\ N — > I + hadrons, where the neutrino disappears and a charged 
lepton I is created, and neutral current (NC) V\ N — > V\ + hadrons, where the neutrino 
survives the interaction. The neutrino-electron scattering cross section is several order 
of magnitude less than neutrino-nucleon cross section and does not affect significantly 
neutrino propagation in matter. The only exception is the W~ resonant scattering for 
the reaction v e e — > W~ + anything, which, however, hardly affects neutrino attenuation 
length in matter for v e energy in the range 2 • 10 15 eV< E v < 2 • 10 16 eV ||. 
The neutrino-nucleon total cross section is approximately <j tot » 10~ 37 cm 2 at E u = 
10 GeV and rises roughly linearly with energy. As a consequence, while in the GeV 
energy range the Earth is completely transparent to neutrinos, it is expected to become 
opaque for suflicently high energy. As to muon neutrinos, the CC interaction lenght inside 
Earth equals Earth's diameter for energy around 40 TeV ||. In order to accurately de- 
scribe the change in energy spectrum of a muon neutrino flux traversing the Earth, both 
the effect of neutrino-nucleon CC and NC must be properly taken into account. The prob- 
lem of neutrino propagation in the Earth can be solved accurately either by Monte-Carlo 
simulation or by approximate iterative solution of the relevant transport equation ||. 
It has been recently pointed out || 0, || |9j that the behaviour of r-neutrinos, whose 
existence should be guaranteed in a neutrino-oscillation scenario, should be significantly 
different from and v e . Whilst muon and electron neutrinos are practically absorbed 
after one CC interaction, the r lepton created by the v T CC scattering may decay in 
flight before losing too much energy, thereby generating a new v T with comparable en- 
ergy. Hence, ultra high energy r-neutrinos should emerge from the Earth instead of being 
absorbed. For a correct evaluation of energy spectrum of z/ T 's emerging from the Earth, 
one has to properly take into account v T interactions as well as r energy loss and decay. 
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An analytical approach similar to that used in Q for has been proposed by S. Iyer et 
al. || neglecting r energy loss. This approach holds as long as energy does not exceed 
10 16 eV because above this energy r interaction lenght becomes comparable with r decay 
lenght (see next section). 

In the present work, a detailed Monte Carlo calculation of v T — r system propagation 
through the Earth has been performed for energy up to 10 20 eV including the r energy 
loss contribution, with special emphasis on the initial v T spectrum deformation as a func- 
tion of the zenith angle of the emerging particle. 



2 Neutrino cross sections and charged lepton inter- 
actions 

The Montecarlo simulation has been performed following neutrinos and charged leptons 
along their path inside the Earth. The Earth model considered is the preliminary Earth 
model of ref. |fT0| . During particle propagation, the occurrence of interactions relevant 
for energy loss, decay or leptons production has been simulated by a unidimensional 
approach. The only secondary particles which have been followed are the u T 's produced 
by t decays and r's arising from v T charged current interactions. 

Deep inelastic neutrino-nucleon scattering is the dominant interaction of energetic 
neutrinos into conventional matter. Charged and neutral current differential cross sections 
used in this work have been calculated in the framework of QCD improved parton model. 
Charged current cross section for a neutrino-isoscalar nucleon interaction reads: 



d 2 a 2G\ME V ( Mfy 



dxdy 7T \Q 2 + Ml 



xq(x,Q 2 )+xq(x,Q 2 )(l-y) 2 , (1) 



where — Q 2 is the invariant momentum transfer between the incident neutrino and out- 
going muon, v = E v — Ei is the energy loss in the lab (target) frame, x = Q 2 /2Mu and 
y = u/E u , M and My/ are the nucleon and intermediate-boson masses, Gp is the Fermi 
constant and the quark distribution functions are: 

q(x,Q 2 ) = Uv{x ' Q2) - dv{x ' Q2) + Us{x ' Q2) ± Q2) + s s (x, Q 2 ) + b s (x, Q 2 ) (2) 

q(x,Q 2 ) = u ^Q 2 )+ ds(X > Q2) +c s (x,Q 2 )+t s (x,Q 2 ), 

where the subscripts v and s label valence and sea contributions, and u, d, c, s, t, b denote 
the distributions for various quark flavors. 

Since a large contribution to cross sections comes from the very low x region laying be- 
yond present accelerator domain, an extrapolation of parton distribution at very low x 
(x > 1CT 8 ) is necessary. We have used the parton distribution set CTEQ3-DIS |TT| 



available in the program library PDFLIB [T2[], with NLO-DGLAP formalism used for Q 



y2 



evolution and assuming for very low x the same functional form measured at x = (9(10 ). 
Even if more sofisticated approaches for low x extrapolation have been developed using 
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dynamical QCD [[L3|] the results for cross section calculations (see Fig. f|) do not differ 
more than 10% from the approach taken here with CTEQ3-DIS plus "brute force" ex- 



trapolation p, 



The electromagnetic interactions of muons with matter, at the actually considered en- 
ergies, are dominated by radiative processes rather than ionisation [|14]]. The cross sec- 
tions for electromagnetic radiative processes of r are lower than muon's, still radiative 
interactions remain the dominant process for r energy loss. The cross sections used 
for radiative electromagnetic interactions of r leptons are based on QED calculation for 



Bremsstrahlung |L5[ , for direct pair production |16| and for photonuclear interactions [|17 



by replacing muon or electron mass with r mass (only formulae with explicit lepton mass 
dependence in the original paper have been considered). For all above processes we have 
implemented stochastic interactions for v = (Ef — E^)j E^ > 10~ 3 and a continuous energy 
loss for v = (Ef — Ei)/Ei < 10~ 3 , where Ei and Ef are the r energies before and after 
the interaction respectively. 

It is worth mentioning that Bremsstrahlung cross section scales as the inverse square 
of lepton mass whereas the direct pair production approximately scales according to 
m e /mi fl8|| . As a consequence, the dominant process of r lepton energy loss is direct 



pair production rather than photon radiation. Along with electromagnetic interactions, 
we also considered weak interactions of r with nucleons although they are almost ne- 
glegible even at E T ~ 10 20 eV. The r decay has been simulated by using the TAUOLA 
package [|T9"||. 



3 Simulation results 

In Figs. 0, |3| scatter plots of initial versus final energy are shown. The initial energy is 
the energy of the neutrino before crossing the Earth while the final energy is the energy 
of neutrino emerging from the Earth. The initial energy distribution has been chosen to 
be proportional to E^ 1 so to have a uniform distribution in logE. Only muon neutrinos 
either interacting with neutral current or not interacting at all are able to cross the Earth 
whereas the rest is absorbed in that the muon arising from a charged current interaction 
has a radiation length much lower than its decay length at any considered energy. 
On the other hand, the r generated from a z/ T charged current interaction loses energy until 
its decay length becomes comparable to its radiation length, thereafter it decays into a 
new v T . Unlike for muons, this occurs around 10 16 — 10 17 eV and, as a consequence, 
CC interacting i/ T 's with E T > 10 17 eV, instead of being absorbed, populate the E T < 
10 17 eV energy region. The result of our Monte Carlo calculation shown in Fig. is 
qualitatively in a fairly good agreement with a similar plot in ref. f7j, though a more 
quantitative comparison could be desirable. The z/ T propagation develops in several v t — t 
transformation steps. The accumulation of z/ T energies below 10 17 eV is modulated by the 
zenith angle of the emerging particle because particle traverse different amounts of matter. 
In Fig. f| a scatter plot of final energy versus 6 zenit h is shown. The input energy spectrum 
has been chosen once again oc E^ 1 while the zenith angle distribution is uniform in order 
to better display of the final energy angular dependence. The final energies cluster around 
0.5 • 10 17 eV for nearly horizontal particles and go down to 10 14 eV for vertically incoming 
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particles. 

We have studied the final energy spectrum of tau and muon neutrinos crossing the Earth 
at three different incident zenith angles for three different initial energy spectra: oc E^ 1 
(see Fig. ||), oc E~ 2 (see Fig. ||) and a possible energy spectrum predicted for topological 
defect decays |2(J (see Fig. |7|). The main feature of z/ r spectrum is indeed the presence of a 
characteristic bump before a steep drop, due to the i/ T regeneration mechanism described 
above. This clearly shows up for hard initial spectra (see Fig. ||D whereas it is hardly 
recognizable for softer initial spectra (see Fig. [7]) and becomes a smooth slope in case 
of steeply falling energy spectra (see Fig. ||). The bump position strongly depends on 
zenith angle because energy loss increases as a function of encountered matter thickness. 
The results obtained in this work are in good agreement with those by S. Iyer et al. ||, 
who took an analytical approach. The main difference between the two results is an 
excess of z/ r 's around the characteristic bump in the ex E~ l energy spectrum in our work. 
This can be easily explained by taking into account the difference in the initial neutrino 
spectrum. We have used a pure oc E -1 spectrum up to 10 20 eV while S. Iyer et al. used 
a modulated cutoff for energies greater than 10 17 eV. Since in a oc E^ 1 spectrum the 
number of neutrinos which are shifted from high energy to the characteristic bump at 
lower energy receives approximately an equal contribution for each energy decade above 
the bump, the agreement is restored. Furthermore, in ref. ||, r energy loss is not taken 
into account. However, the effect of such approximation is almost negligible because in 
the energy region considered therein, below 10 15 eV, r decay lenght is much less than its 
interaction lenght. 



4 Conclusions 

Whilst, as it is well known, the Earth start to be opaque for muon neutrino propagation 
at energies around 10 TeV, a regeneration mechanism relevant to z/ T prevents them from 
being absorbed. This mechanism is ultimately owing to the shorter, with respect to 
muon, r lepton decay length. Nevertheless, the unavoidable radiative r energy loss sets 
an upper bound on the energy of emerging u T 's to 10 17 eV. Therefore, for underwater- 
underice neutrino telescopes, whose sensitivity will cover the whole neutrino spectrum, 
z/ r 's could be the major source of neutrino events above the PeV region. In such detectors 
the charactestics of Fig. || could be also used as one of the signatures of tau neutrinos. For 
atmospheric fluorescence detectors the effective cut-off of 10 17 eV could be below their 
energy threshold. If this is the case, those detectors should be able to detect only almost 
horizontal neutrinos. 
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Figure 1: Total u T CC and NC cross sections on isoscalar nucleon target as a function of 
energy. 
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Figure 2: Scatter plot of initial versus final energy of //^propagating through the Earth 
with an incident zenith angle uniformly distributed in the range 95° < ze nith < 180°. 
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Figure 3: Scatter plot of initial versus final energy of z/ T propagating through the Earth 
with an incident zenith angle uniformly distributed in the range 95° < ze nith < 180°. 
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Figure 4: Scatter plot of final energy as a function of zenith angle for an initial energy 
spectrum oc E~ l . 
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Figure 5: Final energy spectra of i/ T 's an d v ^ crossing the Earth at zeinth angles of 180°, 
140° and 100° for an initial energy spectrum oc E~ x . 
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Figure 6: Final energy spectra of i/ T 's an d crossing the Earth at zeinth angles of 180°, 
140° and 100° for an initial energy spectrum oc E~ 2 . 
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Figure 7: Final energy spectra of u T 's and u^s crossing the Earth at zenith angles of 180°, 
140° and 100° for an initial energy spectrum calculated for topological defect decays in 
ref. [EDI. 
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